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Abstract. We analyse the spectral energy distribution (SED) of U, B, V and i-dropout samples from 
GOODS-MUSIC and we determine their physical properties, such as stellar age and mass, dust attenuation 
and star formation rate (SFR). Furthermore, we examine how the strength of Lya emission can be constrained 
from broad-band SED fits instead of relying in spectroscopy. We use our SED fitting tool including the 
effects of nebular emission and we explore different star formation histories (SFHs). We find that SEDs are 
statistically better fitted with nebular emission and exponentially decreasing star formation. Considering this 
result, stellar mass and star formation rate (SFR) estimations modify the specific SFR (SFR/M*) - redshift 
relation, in comparison to previous studies. Finally, our inferred Lya properties are in good agreement with 
the available spectroscopic observations. 



Keywords: Galaxies: starburst, Galaxies: ISM, Galaxies: high-redshift, Ultraviolet: galaxies 



1 Introduction 



Understanding how galaxies assemble their mass is an important goal in high redshift studies (Stark et al. 
Using new SED modeling techniques taking into account nebular emission QSchaerer fc de Barros||2QQ9 



2009). 



2010), 



we present preliminary results of a study of a large sample of high redshift galaxies. We find that statistically, 
the model which provides the best fit is a model including nebular emission, with an exponentially decreasing 
star formation history. This implies a new estimation on age, stellar mass and SFR. Also, we here demonstrate 
that the properties of Lya can also be inferred from broad-band observations, at least statistically for large 
samples. This allows us, for example, to determine trends of Lya with redshift and other parameters, without 
resort to spectroscopy. 

We adopt a A-CDM cosmological model with Hq=70 km s _1 Mpc -1 , Q m =0.3 and ^a=0.7. 



2 Data and method 



We have used the GOODS-MUSIC catalog of |Santini et ah] ( [2009] ) for the GOODS-South field, providing 
photometry in the J7, .B435, Vqoq^ 1775, ^850, J, K, bands mostly from the VLT and HST, and the 3.6, 4.5, 



5.8 and 8.0 /im bands from the IRAC camera onboard Spitzer. Using standard criteria as in Stark et al. (2009), 



we have selected U, B, V and z-drop galaxies. To reduce the contamination rate, we have only retained the 
objects whose median photometric redshifts agree with the targetted redshift range. We are thus left with a 
sample of 389, 705, 199 and 60 galaxies at z ~ 3, z ~ 4, z ~ 5 and z ~ 6. 



Our SED fitting tool as already being described in Schaerer 
( 2010[ ), here, we briefly summarize. 



de Barros (2009) and Schaerer & de Barros 



We use a recent, modified version of the Hyperz photometric redshift code of Bolzonella et al. (2000), taking 



into account nebular emission (lines and continua) which can impact broad brand photometry and derived 
properties (|Schaerer fc de Barros||2009 , |2010| |Ono et aL 2010). We consider a large set of spectral templates 
( |Bruzual fc Chariot 2003), covering different metallicities and a wide range of star formation (SF) histories 
(exponentially decreasing, constant and rising SF) and we add the effects of nebular emission. 
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We adopt a Salpeter IMF (Salp eter|1955 ) from 0.1 to 100 M , and we properly treat the returned ISM mass 
from stars. IGM is treated following |Madau ( 1995 ) and the extinction is described by the Calzetti law ( |Calzettij 
et al. 1120001). Nebular emission from continuum processes and lines is added to the spectra predicted from the 



GALAXEV models as described in Schaerer & de Barros (2009). We define three models: 



• Reference model: constant star formation, age > 50 Myr and solar metallicity. 

• Decreasing model: exponentially decreasing star formation (SFR oc exp(— t/r)). Metallicity and r are free 
parameters, r varying from 10 Myr to oo (SFR = constant). 



• Rising model: we use the mean rising star-formation history from the simulations of Finlator et al. (2011 ). 
Metallicity is a free parameter. 



Furthermore, we define two options: first, -f NEB which indicates that we include nebular emission, with all 
lines except Lya, since this line may be attenuated by radiation transfer processes inside the galaxy or by the 
intervening intergalactic medium and second, -f NEB+Lya which indicates that we include nebular emission, 
with all lines. In all cases, we consider z e [0,10] in steps of 0.1 and Ay e [0,4] mag in steps of 0.1. 

To determine properly confidence intervals for all the parameters, we ran 1000 Monte Carlo simulations 
for each object by perturbing the input broadband photometry assuming the photometric uncertainties are 
Gaussian. This procedure provides the probability distribution of the physical parameters for each source and 
for the ensemble of sources. 



3 Fit quality 

For the four samples and each SFH (see figure [I] for Rising model, we find similar results for all SFHs), ~35% 
(30%-39%) of the objects are best fitted without taking into account nebular emission. This fraction is inde- 
pendent of parameters like M1500 or the filter's number available. For the three SFH, in J7, V and i-dropout, 
we found respectively 68%, 71%, 71% and 77% common objects best fitted without nebular emission and 85%, 
80%, 94% and 88% common objects best fitted taking into account nebular emission. 

Ha is a strong line at 656.4 nm (reference frame) and must affect 3.5/im-4.6/im color for objects between 



z=3.8 and z=5 ( |Shim et al. 2011). We selected B-dropout objects with 3.6/im and 4.5/im avalaible data 
(excluding non-detections) and with median redshit between 3.8 and 5. We obtain a subsample of 303 objects, 
with similar distribution for the best fits (~35%-~65%). Figure [I] shows that objects best fitted with nebular 
emission are almost all blue objects, which can be easily explain by strong nebular emission. We define two 
categories of objects: Wneb, which are objects best fitted without nebular emission and Sneb, which are objects 
best fitted with nebular emission, choosing arbitrarily the "Decreasing model" to do this selection, considering 
the similarity between these two categories among the SFHs. 

Wneb objects show a slight Xr improvement (15% to 25% lower) for models without nebular emission in 
comparison with the Xr f° r models including nebular emission. At the opposite, Sneb objects show a large 
improvement of the Xr f° r models including nebular emission (36% to 51% lower). In all cases, we find lowest 
Xr with Decreasing/Rising models and Sneb are significantly best fitted with Decreasing+NEB/+NEB+Lya 
models. 



4 Star formation rate and stellar mass 



Considering Wneb and Sneb, if these two populations have respectively intrinsic weak and strong nebular 
emission lines, we expect to find two different SFR estimations (or other difference on physical properties), with 
a larger SFR for Sneb. Actually, we find a difference between these two populations at each redshift. In Figure 
|2j we show for z ~ 4 that Reference model doesn't permit any distinction between the two populations while 
with Decreasing+NEB model, Wneb objects have a median SFR lower (by a factor 4) than Sneb. 

For example, at z ~ 4, considering respectively Decreasing+NEB and Rising+NEB model, median SFR 
is increased by 15% (180%) and the median stellar mass is decreased by more than 50% (60%). A direct 
consequence is that we establish a different sSFR-redshift relation (see Figure [3]) in comparison with previous 
study ( jGonzalez et al.||2010 ). 
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Fig. 1. Left: In blue: best fits distribution among Rising, Rising+NEB+Lya and Rising+NEB model and in green: 
best fits distribution between Rising+NEB+Lya and Rising+NEB model. Right: 3.6/xm-4.5/xm color histogram for a 
sub-sample of z e [3.8, 5] objects. In red, best fitted objects with nebular emission and black, best fitted objects without 
nebular emission. 




Fig. 2. Composite probability distribution of M1500 and SFR for the Reference and Decreasing+NEB model for the 
sample at z ~ 4 as determined for each galaxy from our 1000 Monte Carlo simulations. The points overlaid show the 
median value properties for each object in the sample, black dots for Wneb and white dots for Sneb. The overlaid contour 
indicate the 68% integrated probabilities on the ensemble properties measured from the centroid of the distribution. The 
dashed line represents the Kennicutt relation ( Kennicutt|19 98 ) . 



5 Constraint on Lya 

In some models, we introduce an additionnal free parameter, a variable Lya strength described by the relative 
Lya escape fraction P^ a e [0,1], defined by L(Lya) = f£^ a x L jB (Lya), where L B is the intrinsic Lya luminosity 
of the spectral template given by its Lyman continuum flux and the case B assumption and L(Lya) is the adopted 
Lya luminosity for the spectral template (before any additional attenuation with the Calzetti law, assumed to 
affect stars and gaz in the same manner). Values f£^ a < 1 therefore describe an additional reduction of Lya 
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Fig. 3. Evolution of the specific SFR with redshift. Black squares: Reference model, red dots: Reference+NEB model, 
blue triangles: Decreasing+NEB model and green stars: Rising+NEB model. +NEB models are used for all redshift 
except z no 6, for which we use +NEB+Lya models. 



Overall it turns out that f£^ a is not well constrained for individual objects. For each sample, however, the 
pdf shows two relative maxima close to f£ eZ 



beyond the attenuation suffered by the UV continuum. 

rel 
Lya 

Tya ~ an d I? whose relative importance varies between the samples. 
To quantify this behaviour further, regardless of the detailed shape of the pdf, we count the number of objects 
with EW(Lya) > 50 A. We define the corresponding fraction of Lya objects as iL ya = N(EW(Lya))/Ntot- 
In Figure |4j we plot TLya as a function of the absolute UV magnitude for the sample of z ~ 3, 4 and 5. We 
find that the fraction of objects showing Lya emission increases with redshift and at each redshift, we find 
that Lya enission is more common in galaxies with fainter UV magnitudes. These two main results are in 
good agreement with previous spectroscopic studies ( Stark et al.||2011 Shapley et al.|[2003 Hayes et al.||2011 ), 
showing that photometric data can also provide information on Lya emission. 



6 Conclusions 



Using an updated version of the Hyperz photometric redshift code of Bolzonella et al. (2000) adding nebular 



emission (lines and continua) to the spectral templates (Schaerer & de Barros 2009, 2010), we have analysed 
a large sample of Lyman-break selected galaxies at z ~ 3-6 in the GOODS-S field, for which deep multi-band 
photometry from the [/-band to 8 /am is available. 

We find that at each redshift, an exponentially decreasing star formation history with nebular emission 
provides statistically better SED fits than models with constant or rising star formation. Furthemore, decreasing 
and rising star formation with nebular emission are the only models coherent with 3.6-4.5 /im color which is 
likely influenced by Ha emission at z ~ 4 (Shi m et al.|2011 ). These results have direct consequences on physical 
parameters estimation like stellar mass and SFR, which imply possibly a new mass assembly interpretation (de 
Barros et al. 2011, in preparation). Finally, we show that significant trends of Lya strength with redshift and 
with UV magnitude can be inferred from broad-band photometry observations of large samples of galaxies using 
our models (Schaerer et al. 2011, accepted). 
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Fig. 4. Fraction of galaxies with a large Lya equivalent width EW(Lya) > 50 Acompared to the observed fraction 
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